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ABSTRACT 

Whether protoplanetary disks accrete at observationally significant rates by the magnetorotational 
instabihty (MRI) depends on how weU ionized they are. Disk surface layers ionized by stellar X-rays 
are susceptible to charge neutralization by small condensates, ranging from ^0.01-/im-sized grains 
to angstrom-sized polycyclic aromatic hydrocarbons (PAHs). Ion densities in X-ray-irradiated sur- 
faces are so low that ambipolar diffusion weakens the MRI. Here we show that ionization by stellar 
far-ultraviolet (FUV) radiation enables full-blown MRI turbulence in disk surface layers. Far-UV 
ionization of atomic carbon and sulfur produces a plasma so dense that it is immune to ion recom- 
bination on grains and PAHs. The FUV-ionized layer, of thickness 0.01-0.1 g/cm^, behaves in the 
ideal magnetohydrodynamic limit and can accrete at observationally significant rates at radii > 1-10 
AU. Surface layer accretion driven by FUV ionization can reproduce the trend of increasing accretion 
rate with increasing hole size seen in transitional disks. At radii < 1-10 AU, FUV-ionized surface 
layers cannot sustain the accretion rates generated at larger distance, and unless turbulent mixing of 
plasma can thicken the MRI-active layer, an additional means of transport is needed. In the case of 
transitional disks, it could be provided by planets. 

Subject headings: accretion, accretion disks — instabilities — magnetohydrodynamics (MHD) — pro- 
toplanetary disks — stars: pre-main sequence — ultraviolet: stars 
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1. INTRODUCTION 



T Tauri stars, Her big Ae/Be stars, and yo ung brown 



dwarfs accrete (e.g., Hartmann et al. 2006). Perhaps 



the most direct evidence tor accretion comes from obser 
vations of radiation at blue to ultraviolet wavelengths, 
emitted in excess of the Wien tail of the stellar photo- 
sphere. The ultraviolet excess arises in large part from 
gas that is shock he ated upon free-fallin g onto the stel- 
lar su rface |Calvet fc Gullbring 1998; Johns-KruU et al 



2000 Calvet et al. 12004). From these observations, and 



other s ot optical emission lines (e. g., Hartmann et al 



1994 Herczeg fc Hillenbrand 2008), it is inferred that 



young solar mass s tars accrete at rates M 



10- 



lO^^M© yr^i (e.g., MuzeroUe et al. 2005). Any indi 



vidual syst em may exhibit o rder- unity variations in M 
with time ( [Eisner et aL]|2010 ). 

Strictly speaking, these observations imply only that 
gas in the immediate stellar vicinity — i.e., within a few 
stellar radii, in and around the star's magnetosphere — is 
accreting. Gas and dust orbit at much greater stellocen- 
tric distances, residing in disks that in man y cases extend 

100 AU (e.g 



Watson et al. 



conti nuously from ~0.1 to 

20071, but there is no direct indicati on that any ot this 
material is actually accreting (but see Hughes et al.|2011 
for an attempt at measuring non-Keplerian motions ot 
circumstellar gas). Nevertheless the standard assump- 
tion is that infalling gas near the star was transported 
there through the disk, traveling decades in radius over 

[Electronic address: perez-becker@berkeley.edu| 



the system age of t 



age 10^ yr. This assumption is at 
least consistent with the fact that disk masses at large 
radii — as inferred from continuum dust emission at ra- 
dio wavelengths (e.g. 



Andrews & Williams 2005 1 — are 



similar to Mtago- 

One of the most intensively studied mechanisms for 
transporting material is turbulence driven by t he mag- 
netor otational instability (MRI; for a review, see Balbus 
2009). The MRI requires that gas be ionize d enough 
to cou ple dynamically to disk magnetic fields. |Gammie"l 
( 1996 1 appreciated that for protoplanetary disks, which 



are predominantly neutral, the MRI may not operate 
everywhere. Over a wide range of radii at the disk mid- 
plane, the magnetic Reynolds number 



Re 



cji 
D 



10 



-13 



/ T 



ViooKy VAu/ 



(1) 



is so low that MRI turbulence cannot be sustained 
against Ohmic dissipation. The magnetic Reynolds num- 
ber compares the strength of magnetic induction (flux 
freezing) to Ohmic dissipation. Contemporary numerical 
simulations suggest that the MRI requires Re > Re* « 
10^-10"*, the threshold depending o n the initial field ge- 
ometry (e.g., Fleming et al.| 2000). In equation Q, a 
is the disk radius; 1' is the gas temperature; Cg is the 
gas sound speed; h = Cg/fi is the gas scale height; ^l is 
the Kepler orbital frequency; D = 234 (T/K)!/^ a;-i cm^ 
s^^ is the magnetic diffusivity; and is the fractional 
abundance of electrons by number. 
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Gammie (1996) proposed further that MRI-dead mid- ioni zed surface layers. In a recent attempt to estimate 



plane gas is encased by MRI-active surface layers which 
satisfy the Reynolds number criterion because they are 
ionized sufhciently by interstellar cosmic rays. He es- 
timated that the MRI-active layer extends down to a 
surface density S* « 100 g/cm^ — essen tially the stop 



ping column f or GeV-energy cosmic rays ( Umebayashi & 
Nakano 1981). For comparison, the total surface den- 
sity at a = 1 AU in the minimum-mass solar nebula is 
about 2200 g/cm^ (e.g., [Chiang fc Youdin|[2010| ). Sub- 
sequent work, discussed below, has shown that a fully 
MRI-active surface density of S* 
overestimate. 



Glassgold et al. 



(1997) pointed out that X-rays emit- 
ted by active stellar coronae (e.g., Preibisch et al.|[2005 l 
are more important than cosmic rays in determining the 
ionization fraction, in part because disk surface layers 
are likely shielded from cosmic rays by magnetized stel- 
lar outflows. Such shielding is almost certainly present, 
as its effects are appreciable even for the comparatively 
feeble Solar wind. At cosmic ray energies of 0.1-1 GeV, 
the present-day Solar wind lowers the Galactic cosmic 
ray flux at Earth by roughly an order of magnitude com- 
pared to e stimated flux es outside the heliosphere (see 
Figure 1 of Reedy||1987 |. Moreover, the anti-correlation 



between the terrestrial cosmic ray flux and Solar activity 



has been cle arly observed over many Solar cycles ( Svens- 
mark 1998). By comparison with the Solar wind, pre- 
mam-sequence winds have mass loss rates that are five or- 
ders of magnitude larger, and thus we expect the shield- 
ing of disk surface layers from Galactic cosmic rays to 
be essentially complete for young stellar systems. Nev- 
ertheless cosmic rays might still leak in from the "side," 
entering the disk edge-on from the outside. The impact 
of "sideways" cosmic rays is difficult to assess because 
whether cosmic rays can penetrate radially inward de- 
pends on the radial distribution of matter and on the 
magnetic field distribution, both of which are unknown. 
Cosmic rays might spiral toward the disk along interstel- 
lar field lines only to be mirrored back. If we assume that 
this does not happen, and are willing to be guided by the 
surface density profile of the minimum-mass solar neb- 
ula, we estimate that ionization by sideways cosmic rays 
might be significant on the outskirts of the disk at a > 30 
AU. Under these generous assumptions, cosmic rays can 
sustain MRI turbulence in the entire outer disk — from 
the disk surface to the midplane, at a > 30 AU. We will 
r eturn brie fly to this scena. ri o in Section [4[ 



Bai fc G oodman] p009|, iTurner et al. 



(2007), and 
ayers could 
They found 



Turner et al.,(,2010| ) studied how disk surface 
be made MRI active by stellar X-rays. 
E* ^ 1-30 g/cm^ — the range over which photons of en- 
ergy 1-10 keV are stopped. In their models, the pre- 
cise extent of the MRI-turbulent region depended on the 
abundance of charge-neutralizing grains. X-ray ionized 
layers, though less deep than the layer envisioned by 



Gammie (19961, were considered by Turner et al.| ( |2010 l 
to be capable of transporting mass at observationally sig- 
nificant rates at a = 1 AU. 

All of the aforementioned studies of layered accretion 
assumed that Ohmic dissipation lim its S* , as gauged b y 



Re or the related Elsasser number (Turner et al. 2007) 



This assumption, widespread in the literature, overes- 
timates S* and the vigor of MRI turbulence in X-ray 



E*, Perez-Becker fc Chiang (2011, hereafter PBCll) em- 



ployed an X-ray ionization model similar to those used 
by previous workers, b ut accounted for two a dditional 
factors. First, following Chiang fc Murray-Clay (2007), 
they tested whether ambipoiar dittusion — the"decoupling 
of neutral matter from plasma — limits E*. They gauged 
the importance of ambipoiar diffusion by computing at a 
given depth the number of times a neutral II2 molecule 
collides with charged species in a dynamical time 



100 g/cnr is a gross 



^chargc/^in 



1 



^charge 

10-8 



)( 



ntot 



101 



T-3 



3/2 



(2) 

Here nchargc = sJchargc'^tot IS thc total number density 
of singly charged species (assuming multiply charged 
species are negligible); ^charge is the fractional charge 
density; ntot is the number density of hydrogen nu- 
clei (mostly in the form of neutral H2); and /3in ~ 
2 X 10~^ cm'^ s~i is the coUisional rate coefficient for 
singly ch arged species to share their momen tum with 
neutrals ( [Draine et aLll983| pillar et al.|1997[ ). Whereas 
Re measures how well plasma is tied to magnetic fields. 
Am assesses how well neutral gas is coupled to plasma. 
No matter how large Re may be, unless Am is also large, 
the magnetic stresses felt by the plasma will not be effec- 
tively communicated to the bulk of the d isk gas, which is 
composed overwhelmin gly of neutral H2 ('Blaes fc Balbus 
1994| [Kunz fc Balbus| [2604 : Dcsch 2004; Wardlc ' Mff 
iNumerical simulations ciy jHawley fc StonCj (1998, here- 
after HS) suggested that unless Am > Am* « 10^, MRI 
turbulence would not be sustained in neutral gas. Am- 
bipoiar diffusion is of especial concern in disk surface lay- 
ers because Am depends on the absolute number density 
of charged particles, and not just the relative density as in 
Re. At a fixed ionization fraction, the rapidly decreasing 
gas density with increasing height above the midplane 
renders the MRI increasingly susceptible to ambipoiar 
diffusion. 

The second factor considered by PBCll concerned 
polycyclic aromatic hydrocarbons (PAHs). Like grains, 
these macromolecules act as sites of ion recombination. 
Their abundances in disk surface layers were constrained 
by PBCll using Spitzer satellite observations of PAH 
emission lines, detected in a large fraction of Herbig Ae 
stars surveyed ( [Geers et al 2006 )Q Disk PAHs were 
found to reduce Am and He by factors > 10 (PBCll). 
For typical stellar X-ray spectra. Am reached maximum 
values on the order of unity at E '-^ 1 g/cm^. Such max- 
imum values were attained only for the lowest plausible 
PAH abundances. Comparing max Am ^ 1 with HS's 
determination that Am* ^ 10^, PBCll concluded that 
ambipoiar diffusion, abetted by PAHs, reduced accretion 
rates in X-ray ionized surface layers to values too small 
compared to observations. 

^ Some workers wish to ignore PAHs when studying X-ray-driven 
chemistry on the grounds that they are rarely detected in T Tauri 
stars. It is hard for us to follow this argument. Herbig Ae disks 
commonly evince PAHs and exhibit accretion rates similar to if 
not higher than those of T Tauri disks. Treating the problem of T 
Tauri disk accretion separate from the problern qf^Herbig Ae disk 
accretion seems unjustified. As discussed by ^Geers et ai.j | ^20d6| , 
PAHs usually go undetected in T Tauri stars not because they are 
absent but because they fluoresce less luminously in the weaker 
ultraviolet fields of their host stars. 
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The conclusion of PBCll has since been tempered 
by two new developments. First, PBCll underesti- 
mated Am because they omitted the contribution from 
momentum-coupling collisions between H2 and charged 
PAHs. This neglect is potentially significant because 
charged PAHs are as well coupled to the magnetic field 
as are atomic/molecular ions (HCO+ and metal ions 
like Mg"*" in their simple model). They argued that be- 
cause the abundance of charged PAHs is less than that of 
molecular ions, Am should not be significantly increased. 
This statement applies for the lowest PAH abundance 
considered in PBCll of xp ah ^ 10~^^ (he r e mea sured 



per H), but |Baij (j20TT]) and [Mohanty etall ( |2011[ ) have 
pointed out that PAHs can dominate the charge budget 
at higher Spah- In this high- PAH limit, although PAHs 
decrease the free electron abundance and therefore Re, 
they can boost Am — and in principle M — in comparison 
to the low-PAH case. In this paper we correct for this 
effect by taking 



"charge — ?^i + 'T-PAH(Z=-1) + "PAH(Z=1) 



(3) 



when computing Am, where the various densities refer 
to atomic/molecular ions, PAHs with a single negative 
charge, and PAHs with a single positive charge, respec- 
tively. 



In another recent development, Bai & Stone (2011; 
hereafter BSll) undertook new siniulations of the MRI 
with ambipolar diffusion to update those of HS. They 
confirmed the result of HS that when Am > 10-10^, 
neutrals and ions behave essentially as one fluid with a 
Shakura-Sunyaev transport parameter a that can be as 
high as ~0.1-0.5. Their new result is that MRI turbu- 
lence can still be sustained in the neutrals for Am < 1, 
albeit in a weakened form with max a decreasing with 
decreasing Am. When Am ^ 1(0.1), BSll find that 
max a w 0.01(0.0007). Whether a attains its maximum 
value for a given Am depends on the assumed strength 
and geometry of the background magnetic field. 

In this paper we consider a third, often neglected 
source of ionization: far ultraviolet (FUV = photon en- 
ergies between ~6 and 13.6 eV) radiation emitted by the 
central star. As measured with the Hubble Space Tele- 
scope and FUSE, typical FUV luminosities from young 
stars are > 10^° erg/s, much of it in atomic lines (e.g.. 
T his radiation originates from the 
Gullbring'19981 iJohns 



Bergin et al. |2007[ ). 

stellar accretion shock dCalvet 



Krull et al.||2 000[ |Calv et ct al. 2004 ) and fromThe active 
stellar chromosphere ( [ATexander et ai.||20d5[ ). Photoion- 
ization/photodissociation regions generate d by FUV ra- 



diation have been studied extensively (e.g., Gorti fc Hoi 



lenbachn2008| and references therein; [ l ielens & HolleiT 



Hach 1985), but their implications for the M RI have not 
oeen well quantified. Semenov et al. (2004) carried out 
detailed calculations ot the disk ionization fraction due 
to X-ray and FUV radiation in order to estimate the 
extent of the MRI-active layer, but they did not con- 
sider the charge-adsorbing effects of PAHs or ambipolar 
diffusion. We seek here to give a more comprehensive 
treatment, including new estimates of M driven by FUV 
ionization, made possible in part by the powerful new 
results of BSll. 

High in the disk atmosphere, molecules are photodis- 
sociated and elements take their atomic form. We focus 



on FUV photons (having energies less than the H Lyman 
limit) because they are not stopped by atomic hydrogen 
and as such penetrate more deeply into the disk than do 
Lyman continuum photons. Far-UV photons ionize trace 
atoms. Some elements, when nearly fully ionized, may be 
so abundant that the criteria for MRI accretion are eas- 
ily satisfied. Carbon, for example, has a first ionization 
energy of 11.26 eV and a cosmic number abund ance rel- 
ative to hydrogen of approximately 2.9 x 10"'' (Lodders 
20031 her Table 2). The Stromgren layer of ClI gen- 



erated by FUV radiation is characterized by ionization 
fractions — and, by extension, values of Am and Re — up 
to 10^ times larger than those reported by PBCll for X- 
ray ionized layers of H2. Thus ambipolar diffusion may 
not pose the same threat in FUV-ionized layers that it 
does in X-ray ionized layers. Moreover, electron densi- 
ties in FUV-ionized layers may be so large that they are 
little impacted by PAHs or other small grains. 

Of the various elements that can be ionized by FUV ra- 
diation, we model in this paper only carbon and sulfur. 
Their cosmic abundances are relatively high, and they 
are among the least likely elements to be depleted onto 
grains (as evidenced in the diffu se interstellar medium; 
e.g., Jenkins 2009; Savage & Sembach 1996). Line emis- 
sion from ionizecT carbon (|CII| 158 nm) h as been de- 
tected by Herschel in pr otoplanetary disks ( Pinte et al 
2010 Sturm et al. 2010). The observed line fluxes can 



be reproduced by mod els that assume a co smic gas phase 
abundance of carbon ( [Woitke et al.|[20l0 1 , although our 
own analysis shows that lower abundances are also pos- 
sible because the [CII]-emitting layer is optically thick to 
its own emission. The abundance of sulfur has not been 
so constrained because no lines from sulfur h ave been un- 
ambiguously detected in disks jL ahuis ct al.||2007 Wat- 
'son et al.' 2007; 'Gor ti fc Hollenbach||2008i . INevertheless 
^Mcijerink et al. ( 2008| tound that a sultur abun dance of 
up to 7 X 10~'^ per H, or ~2/5 that of cosmic (Lodders 



2003 1 , could be reconciled with the non-detection of the 
[blj 25 /zm line. The true upper li mit on the sulfur abun 



dance may be even higher because iMeijerink et al. ( 2008 1 
did not account for FUV ionization of SL 

Our paper is structured as follows. The equations gov- 
erning the ionization balance of carbon and sulfur in disk 
surface layers are presented in Section [2j Account is 
made of PAHs; sub-micron sized grains which can attenu- 
ate FUV radiation; and of H2 which can absorb photons 
having energies > 11.2 eV. Results for the carbon and 
sulfur ionization fronts, and the values for Am, Re, and 
S* they imply, are given in Section [3[ The impact of 
the HaU effect (e.g., |Wardle fc Salmeron|[20 ll1) on FUV- 
ionized layers is also considered m Section [3j A summary 
is supplied in Section [ij together with estimates of M as 
a function of disk radius, as driven by FUV radiation, X- 
rays, or cosmic rays. We compare our results to stellar 
accretion rates observed in conventional T Tauri disks 
and transitional disks with inner optically thin holes. 



2. MODEL FOR FUV IONIZATION 
At a given radius a from a star of mass M — 1M( 



0: 



we compute the equilibrium abundances of H2, HI, CI, 
CII, SI, and SII as a function of the vertical column pen- 
etrated by FUV radiation. The key equations governing 
the ratios of H2:H, CLCII, and SLSII are, qualitatively. 
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Rate of dissociation of H2 by FUV photons 

= Rate of formation of H2 on grains (4) 



Rate of ionization of CI by FUV photons 
= Rate of recombination of CII 
with electrons, grains, and PAHs 



(5) 



and 



Rate of ionization of SI by FUV photons 
— Rate of recombination of SII 

with electrons, grains, and PAHs . (6) 

For our standard model we take a — Z AU. This pa- 
rameter and others will be varied to explore their impact 
on the MRI-active column E*. 

2.1. FUV Luminosity 

Two ingredients of our model, II2 and C, are disso- 
ciated and ionized, respectively, by FUV photons hav- 
ing practically the same wavelength range: A ~ 912- 
1109A. A third ingredient, S, is ionized by photons hav- 
ing A < 1198 A. In the combined wavelength interval 
A w 912-1198A we take the stellar luminosity of our stan- 
dard model to be Lpuv = 10^^ erg/s, in continuum pho- 
tons. We take Lpuv to be distributed uniformly over this 
wavelength interval so that approximately (2/3)Lfuv 
can be absorbed by H2, C, and S, while the remaining 
(1/3)I/Fuv can be absorbed by S but not H2 or C. 

Our choice for the continuum luminosity is compat- 
ible with F USE observations of T Tauri stars (Bcrginj 
et al. 2003 see their Figure 1; admittedly their spec- 
tra are extrapolated at A < 950A for TW Hydra and 
A < II50A for BP Tau). Our standard ipuv underes- 
timates the true ionizing luminosity because we neglect 
commonly observed FUV emission lines (but note that 
our wavelength range does not include the powerful H 
Lyman-alpha line at 1216A). It also neglects the large 
contribution fr om the hotter photosphe res of more mas- 
sive stars (e.g., 'Martin- Zai'di et al.|2008 |. We make some 
account for these effects by varyin g Lp uv up to 10^^ erg/s 
in our parameter survey (Section 3.3 1. 



2.2. Total Gas Columns and Densities 

We present our results as a function of the total ver- 
tical hydrogen column density TVtot = A^h + 2A^H2 i mea- 
sured perpendicular to and toward the disk midplane. 
An equivalent measure is the mass surface density E = 
Ntotmii, where mn = 1-7 x 10~^^ g is the mass of the 
hydrogen atom. The local number density of hydrogen 
nuclei is approximated by 



"■tot ~ Ntot/h. 



(7) 



Stellar radiation enters the disk at a grazing angle 
9 ^ 3h/a measure d from the flared disk surface (e.g., 
Chiang et al. 20011. We assume in this work that pho- 
tons travel in straight paths, i.e., we neglect scattering 
of FUV radiation into directions other than that of the 



original beam entering the disk. Thus a beam of radia- 
tion that penetrates a vertical column A^tot has traversed 
a larger column parallel to the incident beam direction of 
and is attenuated according to the latter quan- 
tity, not the former. To order of magnitude, ~ 0.3. 

2.3. Gas Temperature 

Gas temperatures in disk surface layers are set by a 
host of heating and cooling processes. Models in the lit- 
erature account differently for these processes and yield 
temperatures which disagree. At a vertical column of 
A'tot ^ 10^^ cm~^ — roughly where the deepest of the 
fronts we compu te, the SI/SII ionization front, may lie — 
Glassgold et al. (2007, hereafter GNI07) reported a gas 
temperature ot T « 60 K at a « 5-10 AU (see their 
Figure 2) 



By comparison, Gorti 



& HoUenbach 
= 8 AU 



(2008, 
the S ioniza- 



hereafter GH08) found that at a 
tion front occurs at T « 600 K (see their Figures 1 and 
5). One reason for this difference seems to be that GNI07 
omitted photoelectric heating from PAHs, which domi- 
nates thermal balance according to GH08. But the PAH 
abundance assumed by GH08 seems too high, exceeding 
by about an order of magnitude the highest plausible 
abundance inferred by PBCll. Reducing the heating 
rate from PAHs by more than a factor of 10 from its 
value in GH08 (see their Figure 4) would imply that X- 
ray heating dominates thermal balance — as it tends to 
do by default in GNI07. 

In this paper we do not solve the thermal balance equa- 
tions, but instead approximate the disk surface layer 
as vertically isothermal at a temperature drawn from 
GNI07. This is similar to the choice we made in PBCll, 
except that there we were interested in X-ray ionized 
columns A^tot ^ 10^^ cm~^, whereas here we are inter- 

-i22 



< 10^^ cm-^. As 
gas temperatures 



ested in FUV-ionized columns A'tot 
A^tot decreases from 10^^ to 10^^ cm~^ 
in GNI07 rise by factors of ~2-3. At still higher altitudes 
for which A^tot < 10^"'^ cm~^, temperatures rise steeply by 
an order of magnitude or more. As a simple compromise 
we multiply the temperature profile of PBCll by a factor 
of 3 but otherwise keep the same scaling: 



240 



3AU 



-3/7 



K. 



(8) 



Our assumption that the gas is vertically isothermal ig- 
nores the steep temperature gradient at A'tot < 10^^ 
cm^^ where dust-gas cooling is not effective. This ne- 
glect should not be serious as the amount of mass con- 
tained in these uppermost layers is small compared to 
the mass in the MRI-active layer, which we will find is 
concentrated near the C and S ionization fronts, located 
within A'tot ~ 10^^-10^*^ cm~^. Moreover, the locations 
of the fronts are not especially sensitive to temperature, 
as we have verified with a few test cases. 

2.4. Very Small Grains (VSGs) 

Very small grains (VSGs) having sizes of order 0.01 
/xm can attenuate FUV radiation. We parameterize this 
absorption with the variable A'ygQ , the total gas column 
for which VSGs present unit optical depth in the FUVj^ 
We consider A^ySG between 10^^ and 10^** cm~^, taking 

^ Equivalent to the variable cth used by GH08. 
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as our standard value Ny^G ~ ^^^^ '-■^ ^ ■ 
fuse interstellar medium, -/VygQ ~ 10^^ cm~^; its value 
in disk surface layers may be much higher — by 1-3 or- 
ders of magnitude as judged by disk infrared spectral en- 
ergy distributions (see Table 3 of PBCll and references 
therein) — because of grain growth and sedimentation. 

Colliding with ions and faster moving electrons, VSGs 
carry on average a net negative charge. They therefore 
act as sites of ion recombination. To account for such 
recombination, we model VSGs as conducting spheres 
of radius sysG = 0.01 /ini. Their number abundance 
relative to hydrogen is 



a:vsG = "vsc/^tot 



vsg'^'^vsg 



(9) 



In PBCll, we solved for the detailed charge distributions 
of grains using a set of recurrence equations. Solving the 
recurrence equations to obtain self-consistent solutions 
for the electron density, the ion density, and the total 
charge carried by grains/PAHs was necessary because 
in the lightly ionized regions irradiated by X-rays, the 
charge carried by PAHs could be significant compared 
to the amount of charge in free electrons. By compar- 
ison, FUV ionized regions are simpler to treat because 
free electrons are orders of magnitude more abundant 
than grains/PAHs. Thus our analysis here is simplified: 
we assume that every grain has an equilibrium charge 
ZvsG such that it collides with ions as frequently as it 
collides with electrons. That is, .^vSG is such that the 
rate coefficients (units of cm'^ s""'^) aysG.ion and aySG.e — 
as given by equations 9 and 10 of PBCll — are equal 
(see Figures 3 and 4 of PBCll). Under this approxi- 
mation, the rate at which an ion (CII or SII) neutral- 
izes by colliding with VSGs is simply ?^VSGQ!VSGI where 
avSG = avSG.ion = avsG,o- The values of aysG so de- 
rived differ by 50-70% depending on whether the domi- 
nant ion is CII (atomic weight 12) or SII (atomic weight 
32); for simplicity we take a mean value. Given our 
choices for T and s, and assumed sticking coefficients 
between VSGs and electrons/ions of 1, we find that for 
our standard model .^ysG ~ — 1 in units of the electron 
chargej^and the corresponding aysG ~ 2.5x 10~^ cm^/s. 

A few parting comments about our treatment of VSGs: 
first, in deciding how VSGs attenuate FUV radiation, we 
do not make any explicit assumption about the grain size 
distribution. Our proxy for flux attenuation -/Vyg^ is the 
column at which FUV radiation is absorbed by grains of 
all sizes. Flux attenuation by VSGs will play a significant 
role in determining the extent of MRI activity (Section 
3.3 1. Second, in deciding how ions recombine on VSGs, 
we do make an explicit assumption about grain size, 
namely we take sysG — O.Ol/im. Ignoring larger grain 
sizes maximizes the number abundance and geometric 
surface area of VSGs, and thus maximizes their relevance 
for charge balance. With this choice for grain size, ion 
recombination on VSGs is competitive with other recom- 
bination pathways for our standard model (-/VygQ = 10^^ 
cm~^), but is negligible for more dust-depleted models 



VSG 



> 10 



,22 



^ Our Z differs slightly from the true average (Z) calculated by 
PBCll. 



2.5. Polycyclic Aromatic Hydrocarbons (PAHs) 

Ion recombination on PAHs is treated analogously to 
ion recombination on VSGs. Each PAH has radius 6 A; 
an electron (ion) sticking coefficient of 0.1 (1); and thus 
an average charge ^pah ~ —0.05 at a = 3 AU. The 
corresponding ion-PAH rate coefficient apAH ~ 1 x 10^^ 
cm^/s. 

For our standard (fiducial) model, we take the num- 
ber abundance of PAHs relative to hydrogen nuclei to 
be spAH = "■PAH/'T-tot — 10~®. This is their maximum 
plausible abundance, as inferred by PBCll. For this 
PAH abundance, the ion recombination rate on PAHs 
will turn out to be only comparable to the ion recom- 
bination rate with free electrons in FUV-ionized layers 
(the contribution from PAHs was much more significant 
for the deeper and more poorly ionized layers consid- 
ered in the X-ray model of PBCll). Thus we expect the 
thickness of the FUV-irradiated, MRI-active layer to be 
insensit ive t o PAHs, a finding we highlight at the end of 
Section 13.31 

Note that under our simplified treatment of recombi- 
nation on condensates, the net charge of our system is 
not zero: VSGs and PAHs carry a net negative charge, 
which together with the contribution from free electrons 
is not balanced by positive CII and SII ions. How- 
ever the deviation from charge neutrality (which was 
not present in PBCll because that study accounted ex- 
plicitly for collisional charging of condensates by ions 
and electrons) is small. That is, we have verified that 
ZpAnnpAn + ZysCnYSG is much less than the electron 
number density over the domain of our calculation. 

2.6. H2: Photodissociation and Re- Formation on Grains 

Molecular hydrogen is photodissociated by FUV pho- 
tons in two steps: a photon of wavelength 912A < A < 
1109 A sends the molecule into the first-excited electronic 
state, and then in the subsequent radiative decay, there is 
ap = 23% probability that the molecule lands in a v ibra- 



tionally excited state wh ere it becomes unbound (Dal 



garno & Stephens 
dissociation rate o 
let al. (1980) 



1970). Our calculation of the 
H2 rnolecules follows that of 



de JongI 

according to which there are t^l = 60 vibra- 
tionally split line transitions connecting the ground and 
first-excited electronic states. These lines constitute the 
so-called Lyman band. 
The dissociation rate is 



^//(mec) (2/3)L 



FUV 



/3ygG/3H. . (10) 



All lines of frequency z/ (energy hv ~ 12 eV) in the 
Lyman band are assumed to be equally spaced and to 
have equal oscillator strengths / = 4.6 x lO^'^. The 
factor 7re^//(r7icc) is the frequency-integrated cross sec- 
tion of a single line, where e and rUc are the charge and 
mass of an electron, and c is the speed of light. Thus 
Vi^b^s.'^ f / (tooc)]/j^ = tTeff approximates the effective cross 
section of the entire Lyman band, averaged over the in- 
cident broadband spectrum of frequency width The 
factor of 2/3 accounts approximately for the fraction of 
Lpuv (as we have defined it for 912A < A < 1198A) that 
covers the Lyman band (see Section 2.1). 

The factor /3ysG — exp[— iVtot/(fi'A'ygQ)] accounts for 
attenuation of FUV radiation by VSGs (recall that iVtot 
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is the vertical column so that iVtot/^ is the column par- 
allel to the incident beam of radiation). The factor /3h2 
accounts for at tenuation of H2-diss ociating photons by 
intervening H2 ( de Jong et al.|[l980 ): 



JH2 [In(THjv^)] 




erfc 



7rA2 



where erfc is the complementary error function. The core 
of each line is assumed to be Doppler broadened with 
frequency width Ai^d = c^v lc\ then 



(12) 



is a measure of the optical depth at line center along the 
incident beam. The Voigt parameter 



h = 7/(47rAz^D) 



(13) 



compares the natural line width 7 = 1.16 x 10® s ^ to 
the Doppler width, and the non-dimensional parameter 



A = AiyL/(2?yLAt'D) 



(14) 



compares the inter-line spacing to the Doppler width, 
where A^-l = 6 x 10"'^^ Hz is the frequency width of 
the entire Lyman band. Equation (11) applies only for 
TH2 1, a condition which is valid over the domain of 
our calculation. 

Although the photodissociating continuum for H2 and 
the photoionizing continuum for CI nearly overlap in the 
FUV, we can safely neglect shielding of H2 by CI because 
in our model the II2 dissociation front occurs at a higher 
altitude (lower A'tot) than the C ionization front. By the 
same token, shie lding of CI by H2 cannot be neglected 
(see Section 2.7 where it is accounted for by the factor 

^C.Ha)- 

Molecular hydrogen forms when two H atoms combine 
on a grain surfa ce, releasing the heat of formation to 
the grain lattice (Gould & Salpeter 19631. The rate of 
formation is given by 



(15) 



where nu and vy^ are the number density and mean ther- 
mal speed of atomic hydrogen, respectively, a nd 77 ~ 0.2 



is the format ion efficiency on olivine grains (Cazaux & 
Tielens||2010| 



At each depth in our ID model, we solve for nfj/T^H 
using the equilibrium condition (|4[ 



-Rdiss — ^fo 



(16) 



2.7. Carbon: Abundance and Ionization Equilibrium 

We take the abundance of atomic carbon in the gas 
phase to be xc = «c/^tot = 10~'''e, where the dimen- 
sionless factor e accounts for the sequestering of carbon 
into grains. For our standard model, e — 1, which gives a 
gas phase carbon abu ndance similar to that of the diffuse 
interstellar medium (Jenkins 20091. By using [CII] 158 
fim line fluxes from disks as measured by the Herschel 
satellite (Pinte et al. 2010), we estimate a lower bound 
on e of -1/30. 



Equation ([5| for photoionization equilibrium of carbon 
is quantified as 



(2/3)L 



FUV 



'^ClO'bf,c/3c/3c,H2/3vSG 



= ncn (ncarcc,c + ?^PAHapAH + "vsGavsc) , (17) 

where hv w 12 eV is the typical FUV photon energy, 
^ici (ncu = nc — nci) is the neutral (ionized) car- 
bon number density, Uc is the electron number den- 
sity, (Tbt,c = 10~^^ cm^ is the b ound-free cross see 



tion for CI (Cruddace et al. 19741, and aroc,c — 5.3 x 
10~^2(r/240 K)~'^-'' cm ^s is the radiativ e recombination 
rate coefhcient for CII ( Pequignot et al.||l991 ) . The fac- 
tor of 2/3 accounts f or th e fraction of Lftiv (as we have 
defined it in Section 2.1) covered by the photoionizing 
continuum for CI. The factor /3c = exp(— A^ciCbf.c/^) ac- 
counts for self-shielding of carbon, where iVci is the verti- 
cal column of CI. The factor /3c, H2 accounts for shielding 
of carbon by the Lymaii band transitions of molecular 
hydrogen ( de Jong et al.|[T980 ): 



1 



7rA2 



exp 



TH2& 
"ttA^ 



(18) 



This factor drops significantly when the H2 column 
N'n.JO > 10^2 cm~2 — large enough for the Lorentzian 
wings of each Lyman line to blacken the entire Lyman 
band. As noted previously, all the FUV photons that 
can ionize CI (1104A > A > 912A) can also dissociate 
H2 (IIO9A > A > 912I). This is true for C but not for 
other elements like S. 

We neglect shielding of CI by SI because the S ion- 
ization front lies at a greater depth than that of the C 
ionization front (because S is a factor of 10 less abundant 
than C everywhere). 

2.8. Sulfur: Abundance and Ionization Equilibrium 

The abundance of gas-phase atomic sulfur is xs = 
ng/ntot = 10~^e. For our standard model, e — 1, which 
yields a sulfur abundance similar to that of the diffuse 
interstellar medium (Jenkins 2009). For simplicity, we 
assume that the depletion factor e for sulfur is the same 
as that for carbon, so that xs/xc = 0.1 for all our mod- 
els. The lower abundance of sulfur relative to carbon im- 
plies that S-ionizing radiation will penetrate to greater 
depths than C-ionizing radiation (see, e.g., the order-of- 
magnitude estimate in PBCll), implying a deeper and 
more massive MRI-active layer. 

Equation ^ for photoionization equilibrium of sulfur 
is quantified as 



(1/3)L 



FUV 



"-sifbt,s/3s/3vsG 



= "SII ('^Garoc,S + "PAHttPAH + "•VSG^VSg) , (19) 



cm^ is the 



where hi' 11 eV, Cbf.s ~ 5. 5 x 10 

photoionizing cross-section for S ( Cruddace et al.|[l974 ) 

a,^^ g = 4.9 X io-i2(r/240K)-° fa^ cm^/ s is the radia- 



tive recombinati on rate coefhcient for SII ( Aldrovandi & 
Pequignot||1973| [Gou ld 1978), and ngi (ngii = ng 



_ nsi) 
is the neutral (ionized) sulfur number density. The fac- 
tor /3g = cxp(— A^gicTbf.s/^) accounts for self-shielding of 
sulfur, where iVgi is the vertical column of SI. 
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The factor of 1 /3 accounts for the fraction of Lpu y (as 



2.11 



we have defined it for 912A < A < 1198A in Section 
that can ionize S but not be absorbed by CI and H2 
TechnicaUy there should be a second term on the left- 
hand side of (19) that is proportional to (2/3)iFuvi to 
account for radiation that can ionize S and be absorbed 
by CI and II2. But this term can be dropped because 
it does not significantly influence the location of the S 
ionization front; the term is negligibly small there, at 
a greater depth than the C ionization front. Note that 
the governing equation for sulfur is coupled to the other 
equations for C and H2 only via the electron density tIq = 
ncn + nsu. 

2.9. Numerical Method of Solution 

Our calculation starts at A^tot = 10^® cm~^, under 
the assumption that the incident flux L-puy/Aira^ is not 
attenuated at this column. We take 10^ logarithmi- 
cally spaced steps to TVtot = 10^"^ cm~^. At each step, 
we solve for the ratio nn/nna using the balance con- 
dition ( jlG] ). The self-shielding factor /3h2 depends on 
ano^so we keep a running integral of uh^, i.e., 

= Z]i("-H2Mot)i(AA'"tot)i- 

By approximating rie = max(?T,cii, ngn), we may solve 
analytically for the ionization states of C and S. At low 
A^tot , nearly all C and S are ionized and ric — ?^CII with 
an error of xc/xg ~ 10%. Equation (17 1 becomes a 
quadratic for ricii/^cij ^nd ?t-.stt /"^ st is determined by 
substituting ncii for rie in equation (19 1. 

At iVtot past the C ionization front, SII becomes the 
dominant ion and we set = nsu- Now equation (19 1 
becomes a quadratic for nsn/ngi. The resultant value 



for ngn is substituted for in equation (17) to solve for 
ncn/nci- 

At each A'tot, we keep track of the running 
columns iVci = X]i("ci/'^tot)i(AiVtot)i and TVsi = 
X]i("si/'^tot)i(AA^tot)i to compute the self-shielding fac- 
tors /3c and /3s, respectively. 

3. RESULTS FOR MRI-ACTIVE SURFACE DENSITIES 



In Section [3T] we describe how the Il2 :H, CLCII, and 
SLSII ratios vary with depth. In Section [X2] we estimate 
the surface densit y S * of the layer that can be MRI- 
In Section [373 we perform a parameter study of 

The MRI-active 



active, 
how S 



varies wit 



1 system properties, 
surface densities E* reported in Section [3T3| are those re- 
sulting from FUV ionization only. In Section [374| we check 
that the Hall effect is not significant for FUV-ionized sur- 
face layers. All our calculations neglect turbulent mix- 
ing of plasma from disk surface layers into the disk inte- 
rior. In Section [33) we try to assess whether this neglect 
is justified, by comparing the dynamical timescale with 
the timescale over which the FUV-irradiated layer equi- 
librates chemically. 

3.1. Photodissociation and Ionization Fronts 

In the top panel of Figure[lJ we show the relative abun- 
dances of species as a function of vertical mass column 
E. High in the atmosphere, hydrogen is mostly in atomic 
form and nearly all of the carbon and sulfur are ionized. 
Deeper down, starting at S ~ 3 x 10~^ g cm~^, atomic H 
gives way to H2. At the H2:H front the Doppler broad- 
ened cores of the various H2 Lyman transitions have be- 
come optically thick. Far-UV photons in the wings of the 
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Fig. 1. — H2 photodissociation front, and C and S ionization 
fronts, for our standard model at a = 3 AU. Results are given as 
functions of tlic vertical surface density S, measured perpendicular 
to the disk midplane (the column actually penetrated by radiation 
parallel to the incident beam direction is T,/9). The correspond- 
ing number density of H nuclei is given on the top x-axis. Upper 
panel: Normalized abundances of atomic H, H2, CII, CI, SII, and 
SI. The H2 photodissociation front occurs at the highest altitude 
above the midplane, and the S ionization front occurs at the lowest. 
Middle panel: Optical depths to FUV radiation presented by H2 
(averaged over the Lyman band), CI, SI, and dust (VSGs). Lower 
panel: Shielding factors by which FUV radiation is attenuated: 
self-shielding of H2 {Pn^); self-shielding of C {Pc)\ shielding of C 
by H2 (/3c, H2); self-shielding of S (/3s); and shielding by very small 
grains (/3vsg)- 
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Fig. 2. — Magnetic Reynolds number Re and ambipolar diffusion 
number Am versus vertical surface density S for standard model 
parameters and a = 3 AU. 

Lyman transitions can stream past the H2:H front and 
continue to ionize C and S at greater depths. 

The CLCII ionization front is located at S ~ 3 x 10"'^ 
g cm~^, where the optical depth to carbon ionizing pho- 
tons is of order unity. At this front, both shielding of car- 
bon by H2 (quantified by /3c. H2) carbon self-shielding 
(/3c) are becoming significant; see the bottom panel of 
Figure [T] 

Far-UV photons having 1198 A > A > 11091 inter- 
act with neither H2 nor C, but can ionize S. These pho- 
tons penetrate both H2 and C fronts, and ionize S at 
greater depths. The SLSII ionization front is located at 
E ^ 1 X 10~^ g cm^^. In this fiducial model, the loca- 
tion of the S ionization front is determined more by flux 
attenuation by very small grains than by self-shielding 
by sulfur; at the front, /3vSG ~ 0.1 and /3s ~ 1. Lower 
but still observationa lly realistic grain abundances are 
considered in Section 3.3 where we will find that the S 




10 



E [gem ^] 



ionization front can be as deep as E ~ 1 x 10 g cm 



Fig. 3. — Maximum transport parameter a and mass accretion 
rate M versus vertical surface density S for standard model pa- 
rameters and a = 3 AU . The value of maxo is computed using the 
empirical relation l |20| of BSll. A vertical dotted line marks the 
MRI-active surface density E* at which, by definition, M peaks. 
The value of S* in FUV-ionized layers is essentially set by the SLSII 
ionization front. It can be as large as ~0.1 g/cm^ for parameter 
choices other than those assumed here; see Figure |4] 

3.2. MRI-Active Surface Density E* 

We measure the extent of the MRLactive column by 
means of the magnetic Reynolds number Re and the 
ion-neutral collision rate Am. Figure [2] shows both di- 
mensionless numbers as a function of E for our standard 
model at a = 3 AU. 

Above the CLCII ionization front, at E < 3 x lO^'^ g 
cm^^, the value of Re remains constant at ~10^° because 
Xc saturates at 10"'* where nearly all the carbon is singly 
ionized. In this region. Am oc E because Am cx ntot- The 
value of Am peaks at ^3x 10'^ near the C ionization front. 
Both Am and Re fall just past the C ionization front, 
but not by more than an order of magnitude: S, which 
is 10 X less abundant than C, remains fully ionized out 
to the S ionization front at E ~ 10~^ g/cm^. Past the S 
ionization front, virtually all FUV photons are used up; 
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the charge density, and by extension both dimensionless 
numbers, drop sharply. 

We define the MRI- active surface density E* as follows. 
Bai & Stone (2011) find from an extensive series of nu- 
merical simulations that for a given Am, the maximum 
value of a is given by 



max a 



50 



-1/2 



(20) 



From the value of max a so computed, we evaluate the 
mass accretion rate 

kT 

A'f = 2 X SttEv = 67rS x max a x — (21) 

as a function of E, where v = acgh is the viscosity, 
^ — Ax 10^^** g is the mean molecular weight, k is Boltz- 
mann's constant, and the prefactor of 2 accounts for the 
top and bottom surfaces of the disk. Note that max a is 
a function of Am, which in turn is a function of E that 
we have computed (Figure [2]). We define E* as that value 
of E for which M peaks; see Figure [sj which plots equa- 
tion (21). According to this definition, E* « 1 x 10^^ 
g/cm . At this column, according to Figure [2] Re re- 
mains larger than Re* , demonstrating that ambipolar 
diffusion limits MRI activity more than Ohmic dissipa- 
tion does, although the effect is slight (the dominance 
of ambipolar diffusion over Ohmic dissipation was much 
more pronounced for the X-ray ionized layers analyzed 
by PBCll). Referring back to Figure [Tj we see that E* 
corresponds essentially to the SLSII ionization front. 

3.3. Parameter Survey, Including Sensitivity to PAHs 

In Figure |4] we explore the sensitivity of E* to various 
parameters. For a given solid curve in a given panel, only 
the one parameter on the x-axis is varied while all others 
are held fixed. Labels on a curve indicate the values to 
which certain parameters are fixed when computing that 
curve. If a parameter is not labeled, its value was set to 
the standard value. 

Over most of parameter space, the MRI-active surface 
density E* in FUV- ionized layers is of order 10~^-10~^ 
g/cm^. The MRI-active surface density is set by the 
depth of the sulfur ionization front. Very small grains 
(VSGs) reduce the depth of this front by absorbing FUV 
radiation (and also by enhancing the ion recombination 
rate, although this is not a dominant effect for FUV- 
ionized layers where the electron density and thus the 
radiative recombination rate are high). At their stan- 
dard abundance — which is probably near their maximum 
abundance (Section |2.4[ ) — VSGs shield sulfur from ioniz- 
ing radiation more than sulfur shields itself; compare /3s 
and PvsG in Figure [l] 

As Ny'^Q increases from 10^^ to 10^"* cm^^ (i.e., as 
the grain abundance decreases), E* increases by a fac- 
tor of 3 (Figure lib). Now most of the shielding of sul- 
fur from FUV radiation is provided by sulfur itself, not 
by VSGs. When grains are insignificant (iVygQ = 10^"* 
cm~^), E* depends on Lpuv a-nd the carbon/sulfur abun- 
dance e more strongly than in the case when grains are 
significant absorbers of FUV radiation (Figure [4];, |4[i). 
The dependences in the case of low grain abundance re- 



semble those calculated by PBCll in their simple FUV- 
Stromgren model: E* cx L^^y and E* cx 1/e. 

The MRI-active thickness E* in FUV-ionized layers 
hardly varies with disk radius a (Fig;ure|4^), because the 
decrease in ionizing flux with increasing a is compensated 
by the lengthening dynamical time (to which Am, 
the factor determining E*, is proportional). 

Finally, E* is not sensitive to the abundance of PAHs 
(Figure pe). For the large electron fractions Xc ~ 10~^- 
10^"* generated by FUV ionization of carbon and sulfur, 
the primary channel for charge neutralization is radiative 
recombination with free electrons. Only at the highest 
PAH abundances does the ion recombination rate with 
PAHs (xpah^pah) approach the ion recombination rate 
with electrons (xcaioc,s)- By contrast. X-ray ionized gas 
has much lower electron fractions and is much more sus- 
ceptible to charge recombination on small condensates 
like PAHs (PBCll). 

The last point is echoed in Figure [5] which shows how 
the free electron fraction varies with PAH abundance in 
both FUV-ionized and X-ray-ionized surface layers. Each 
curve is computed at a fixed surface density characteristic 
of each layer: E « 0.01 g/cm^ for the FUV- irradiated 
layer, and E pa 0.3 g/cm^ for the X-ray-irradiated layer. 
For the FUV-ionized layer, a;o remains constant over the 
range of plausible PAH abundances lO"^'^ < a^PAH ^ 
10^^. For the X-ray- ionized layer, Xc varies by about 
two orders of magnitude over the same range of PAH 
abundance; Xc oc I/ccpah for xpah ?J 5 x 10~^'^ (see 
Section 3.2.1 of PBCll). 

3.4. Hall Diffusion 
When comp uting S*, we have ignored the effects of 



Hall diffusion jWardle| njOO) |Balbus fc Terquem| [2001 
Sano & Stone [20U2[ ). When Hall diffusion dominates 



only electrons are coupled to magnetic fields, and ions 
are de-coupled from magnetic fields by collisions with 
neutrals. Under these circumstances, the evolution of 
the MRI depends on the direction of the magnetic field B 
with respect to the angular frequency O. Hall diffusion 
can increase/decrease E* by an order of ma gnitude or 
more when B a nd H are parallel/anti-parallel (Wardle & 
Salmeron||2011[ ). 

'lb assess the importance of Hall diffusion in FUV- 
ionized surface layers, we evaluate the Hall parameter 



Ha = 



eBxi 
rjHa^ mucVl 



(22) 



Here xi = Jij/ntot is the fractional abundance of ions 
of density n; relative to hydrogen nuclei of density 
fHot \ — B / ^47rntot"iH is the Alfven velocity; e 
is the electron charge; c is the speed of lig ht; and 
VHa = cB/Ojiren i) is the Hall diffusivity (e.g., |Wardle 
fc Salmeron||20li] but note that their Hall parameter is 
the inverse of ours, and they assume rij = rie). The di- 
mensionless Hall parameter is the ratio of the inductive 
term to the Hall term in the magnetic induction equa- 
tion. If Ha ^ 1, then Hall diffusion is not important. 

To evaluate Ha, we estimate B from M by assuming 
that the inequality 



B^ > Bl + Bl> 2BrB,, 



(23) 
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Fig. 4. — How S* varies with (a) stcllocentric distance a; (b) our proxy A^ygQ for the abundance of very small grains (see equation 

t; (c) FUV luminosity Lpuv; (d) the gas-phase abundance of carbon and sulfur e (normalized so that e = 1 corresponds to near-solar 
undances); and (e) the abundance of charge-absorbing PAHs. For a given solid curve in a given panel, only the one parameter on the 
X-axis is varied while all other parameters, unless otherwise annotated on the figure, are held fixed at their standard values. The upper 
bounding curves in panels (a), (c), (d), and (e) are calculated assuming Ny"^^ = 10 cm~^, the lowest grain abundance considered. In 
panel (a), the upper bounding curve ass umes e = 1/30. which we estimate to be the lowest possible value still consistent with observations 
of far-infrared [CII] emission from disks | |Pinte et al.|20l"o| [Sturm et al.|2010[ |. A dot in each panel marks the standard (fiducial) model. 



saturates, and that the Maxwell stress which drives ac- 
cretion is given by 



~2h 



(24) 



(see, e.g., Bai fc Goodnian|2009 the factor of 2 in equa- 



tion arises because the disk has a top and bottom 
face). 

In Figure |6] we show that Ha ^ 1 at E* for all a, over 
the entire parameter space that we have explored. Thus 
we conclude that Hall diffusion is not a concern in FUV- 
ionized surface layers (the same may not be t rue for the 
more poorly io nized X-ray-irradiated layers; Wardle ^ 
Salmeron||2011[ ). 



3.5. The Possibility of Turbulent Mixing: Chemical 
Equilibration Timescale vs. Dynamical Timescale 

Throughout this paper, we have computed ionization 
fractions in a static atmosphere. But the MRI-active 
layer is not static; it is turbulent. We might have under- 
estimated E* because turbulence can mix plasma ver- 
tically toward the midplane, deeper into the disk in- 
terior dlnutsuka fc Sano | |2005| |ll gner & Nelsonl 12006 
'I'urner et ai. 20U7[ ) . Mixmg would be ettective it the 



vertical mixing time is short compared to the timescale 



over which ionized layers achieve chemical equilibrium. 
Without a full-out simulation of MRI turbulence, we 
approximate the vertical mixing time as the dynami- 
cal time tdyn — ■ For the chemical equilibration 
time, we substitute the radiative recombination time 
iroc = l/('^oarcc,s)- The ratio irec/idyn is plotted in Fig- 
ure [7](cf. Figure 11 of PBCll). Because W/idyn > 1 
over some portion of parameter space, mixing might well 
be significant in FUV-ionized layers, and might increase 
Yi* above the values we have computed in this paper. 
Whether the increase would be large or small is hard to 
say. On the one hand, iroc/idyn is never far above unity, 
suggesting that turbulent mixing will merely introduce 
order-unity corrections to our estimates. On the other 
hand, as turbulent mixing dilutes the electron density 
TT-c, the timescale ratio ircc/^dyn might increase with in- 
creasing depth, and the process might run away. 

4. SUMMARY AND IMPLICATIONS FOR DISK 
ACCRETION 

Circumstellar disk material must be sufficiently ion- 
ized if it is to accrete by the magnetorotational insta- 
bility (MRI). We have considered in this paper ioniza- 
tion by stellar far-ultraviolet (FUV) radiation. Although 
FUV radiation cannot penetrate the disk as deeply as 
can X-rays, it generates ionization fractions orders of 
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Fig. 5. — Electron fraction as a function of PAH abundance for 
FUV and X-ray ionized surface layers. The FUV curve is computed 

Lpuv = 10^° erg/s, 
; cm"-^. The large 

abundance of electrons generated in FUV-ionized layers is immune 
to the effects of charge recombination on PAHs, over the range of 
plausible PAH abundances (shaded in grey; PBCll). We contrast 
this behavior with the X-ray curve taken from Figure 7 of PBCll 
(a = 3 AU, Lx = 10^^ erg/s, metal abundance xm = 10"*, S = 0.3 
g cm"'^), which shows that PAHs can reduce electron fractions in 
X-ray-irradiated layers by two orders of magnitude. 



for our standard model parameters (a = 3 AU, 
10^2 cm-2, e = 1) at S 0.01 



,VSG, 




Fig. 6. — Hall parameter Ha as a function of a in FUV-ion ized 
surf ace layers. Values of B are computing using equations (|23|l 
and ( |24[ l. The upper bounding curve is computed for our standard 
(fiducial) parameters, at S* ^ 0.01 g/cm^. The lower bounding 
curve is computed for a model depleted in dust, carbon, and sulfur, 
for which S* 0.2 g/cm'^ (see also Figure [4|i). Hall diffusion is 
not significant for FUV-irradiated surface layers. 

magnitude larger. Carbon and sulfur may be the princi- 
pal sources of free electrons and ions, as these elements 
are cosmically abundant and least likely to be depleted 
onto grains. Far-infrared searches for carbon and sulfur 
emission from disks are so far consistent with gas phase 
abundances for both elements within a factor of ~2 of 
solar. Because the electron and ion abundances gener- 




FlG. 7. — Ratio of the ion-electron recombination timescale t^cc = 
(neOrec.s)"^ to the dynamical time tdyn = as a function of a. 

The lower bounding curve is computed for our standard (fiducial) 
parameters, at S* 0.01 g/cm^. The upper bounding curve is 
computed for a dust-depleted and carbon/sulfur-depleted model, 
for which S* a: 0.2 g/cm^ (see also Figure|4|. 

ated from ionized carbon and sulfur are so large, ion- 
ization fractions in FUV-irradiated layers are little im- 
pacted by PAHs. In FUV-ionized layers, fractional ion 
abundances Xi ~ 10~^-10~'*, dwarfing PAH abundances 
of xpAH ~ 10~^^~10^®. More to the point, ion recombi- 
nation on PAHs in FUV-ionized layers is at most com- 
petitive with ion recombination with free electrons. This 
immunity to PAHs does not apply to X-ray ionized lay- 
ers where Xi < 10~^ and where PAHs o r very small (0.01 



2009 Perez-Becker fc Chiang 2011 PBCll). it is re 



lim sized) grains can suppress the MR I (iBai & Goodman 

■ ■ ^BCT 



freshing that l-'UV ionization is robust against the usual 
difficulties plaguing other, weaker sources of ionization. 

In Figure |8] we compute the possible ranges of disk ac- 
cretion rate M due to either FUV or X-ray ionization, 
and compare them against the range of observed stellar 
accretion rates. The accretion rate M driven by FUV 
ionization is derived using equations ( pO| ) and ( |2T| ) for 
S = S* (taken from Figure l4k) and the corresponding 
value of max a. The accretion rate M driven by X-ray 
ionization is computed similarly, using the abundances of 
charged species of the standard model of PBCll, their 
temperature of T = 80(a/3AU)~^/'^ K, and Equations 
([2])-([3]) of this work to compute Am{T,). Although our 

estimates of M driven by X-ray ion ization may have to 
be rev ised because of the Hall effect ( Wardle &: Salmeron 
201l| ), our estimates of M driven by FUV ionization 
should not be: FUV-ionized, MRI-active layers behave 
in the ideal magnetohydrodynamic limit. 

Modulo the impact of the Hall effect on X-ray-driven 
MRI, from Figure[8]we conclude that accretion rates from 
FUV-ionized surface layers are of the same order of mag- 
nitude as accretion rates from X-ray ionized layers. At 
large radii a > 10 AU, FUV-driven accretion rates may 
exceed X-ray-driven rates. At a < 10 AU, their contri- 
butions may be more nearly equal. At a > 10 AU, the 
FUV-irradiated surface layer can sustain accretion rates 
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Fig. 8. — Accretion rates for conventional (hole-less) disks. The 
shaded region labeled "FUV" (bounded by heavy solid lines) corre- 
sponds to the range of accretion rates possibly driven by the MRI 
in FUV-ionized surface layers. The range of FUV-driven accretion 
rates corresponds to the range of MRI-active surface densities S* 
shown in Figure |4k, which in turn corresponds to a range of possi- 
ble abundances for grains, carbon, and sulfur (PAHs, whose abun- 
dances for all our FUV models are set to their maximum value, are 
not significant for FUV-ionized layers). We compute X-ray-driven 
accretion rates using our corrected definition for Am (Equations 
W^, for the low-PAH and high-PAH cases considered by PBCll. 
T^ne case with no PAHs is shown for comparison only. Note how 
little the X-ray-driven M varies, despite the PAH abundance hav- 
ing changed by three orders of magnitude; the small variation 
in M follows from Equation tell, in which reductions in rt; from 
charge-neutralizing PAHs are offset by the increased abundances 
of charged PAHs themselves. In fact, the high-P AH M even ex- 
ceeds slightly the low-PAH M (see also Bai|2011 1. At a < 10 AU, 
accretion rates in FUV-ionized layers may be cornparable to those 
in X-ray ionized layers. At a > 10 AU, FUV-driven accretion rates 
tend to be larger than X-ray-driven rates. The light shaded region 
labeled "Observed Rates" brackets the range of stellar accretion 
rates for stars of mass hi 0.3-1Mq as shown in Figure 5 of 
[Muzerolle et al.| | )2005j . 

similar to those observed. At a < 10 AU, surface layer 
accretion rates, driven either by FUV or X-ray ionization, 
can still be observationally significant, but they diminish 
with decreasing radius. At a < 1 AU, surface layer ac- 
cretion rates fall below the range typically observed for 
young solar-type stars. 

The problem of too low an accretion rate at small ra- 
dius might be alleviated by turbulent mixing of plasma 
from FUV-irradiated disk surface layers into the disk in- 
terior, as such mixing would enhance the thickness of the 
MRI-active layer. Our crude estimate of the timescales 



involved (Section 3.5 1 suggests that this possibility is 
worth further consideration. 

Note in Figure 8 how the X-ray-driven M actually in- 
creases from the low-PAH case to the high-PAH case. 
This surprising effect arises because increasing the PAH 
abundance (over this particular range from xpah — 
10^^^ to 10~®) increases the abundances of charged 
PAHs themselves, which more than offsets the decreased 
abundances of free electrons and ions. In fact, in the 
high-PAH case, positively and negatively charged PAHs 
are approximately equal in number and represent by far 
the most abundant charged species. The primary recom- 




a [AU] 



Fig. 9. — Same as Figure [s] but showing disk accretion rates 
driven by the combined effects of X-ray and cosmic ray ionization. 
Galactic cosmic rays cannot reach the disk through its top and 
bottom faces because of shielding by the magnetized stellar wind. 
They can, however, penetrate the disk "sideways" through the out- 
ermost portions of the disk, parallel to the midplane — assuming 
they are not magnetically mirrored away. To generate this plot 
we assume a total (vertical) surface density of 2200(a/l AU)~^'^^ 
g/cm^, and a cosmic r ay ionization rate Co exp(— Sg/Sp), where 
Qq = 1/4 X 10"" s-i ( |CaseUi et al.|l99 8!) and Sn = 96 g/cm^ 



the stopping column for Lie V -energy cosmic rays ( [Umebayashi &[ 
|Nakano|198l|l . The factor of 1/4 in Co is our estimate tor the trac- 
tion ot tne celestial sphere (centered at the midplane) not shielded 
by stellar winds. The attenuating column So is obtained by ra- 
dially integrating the volume mass density at the midplane from 
a to infinity. For a > 30 AU, sideways cosmic rays may provide 
enough ionization to render the entire disk MRI active (PBCll). 
Sideways cosmic rays cannot penetrate the disk at a < 10 AU be- 
cause of intervening disk material, and so the accretion rates there 
are practically identical to those for the X-ray-only case (plotted in 
Figure^. The solid and dotted curves correspond to the cases of 
low andnigh PAH abundances, respectively. The main assumption 
underlying this plot is that sideways cosmic rays are not magneti- 
cally mirrored away; this is an uncertain prospect. 

bination pathway is positively charged PAHs colliding 
with negatively charged PAHs (reaction 14 in PBCll). 
Thus as the overall PAH abundance increases, charged 
PAHs can increase Am at a given column E. This effec t 
was missed by PBCll, and is highlighted by Bai (2011 ). 



The consequent enhancement in M is mitigated by the 
decrease in the free electron fraction and thus in Re. In 
computing the X-ray-driven M for the high-PAH case 
in Figure 8, we accounted for the mitigating effects of a 
lower Re by evaluating Am at the E for which Re « 100, 
below which Ohmic dissipation would weaken the MRI. 

One concern is whether gas pressures in FUV-ionized 
layers are lower than the magnetic field pressures re- 
quired to drive our computed accretion rates. If the 
pressure ratio 



Y.*kT/{^ih) 



(25) 



is < 1, magnetic tension defeats the MRI. An upper 
bound on /3piasma is obtained by substituting the lower 
bound on inferred from equations (23 1 and (24). For 
a representative case E* — 0.1 g/cm"^, equations (23)- 



( 25 ) combine to yield the topmost slanted line in Figure 
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Fig. 10. — Same as Figurejslbut for transitional disks and show- 
ing only ou r FUV model. Observed rates are taken from Cal vet] 
|et al.|(|2005"l diamond s), jEspaillat et al.,t2007,.2008, 2010. squares77 
and |KTm et al.| l |2009[ circles with error barsj . I'h eore tical rates for 
FUV -ionized layers are computed with equation \26\ , with S* re- 
computed using 9 ~ 1. The observed trend of increasing Mtrans 
with increa sing aj-im is reproduced (see, however. Section 4.2 of 
IKim et al.|[2b09i for how the observed trend could reflect other 
correlations between stellar mass M and a^im, and M and M). 

[8] above which max /3piasma < 1 and the MRI cannot op- 
erate. The range of FUV-driven accretion rates we have 
computed sits safely below this line. 

In sum, surface layer accretion by FUV ionization can, 
by itself, solve the problem of protoplanetary disk ac- 
cretion at large radius, but not at small radius (unless 
turbulent mixing of plasma can substantially thicken the 
MRI-active layer) . This statement remains one of princi- 
ple and not of fact, because MRI accretion rates depend 
on the strength and geometry of t he background mag- 
netic field threadi ng the disk (e.g., Fleming et al.|[200Ci{ 
Pessah et al.|2007 Bai & Stone 2011], and unfortunately 



the field parameters are not known for actual disksj^ 

For completeness, we show in Figure [9] how sideways 
cosmic rays can, in principle, enhance accretion rates in 
the outermost portions of disks (see the Introduction). 
For an assumed radial surface density profile resembling 
that of the minimum-mass solar nebula, accretion rates 
at a > 30 AU can be strongly increased by sideways cos- 
mic rays. Although cosmic-ray induced accretion rates 
at large radius can be competitive with FUV-induced 
rates, it is not clear that cosmic rays are not mirrored 
away from disks by ambient magnetic fields. Far-UV 
ionization does not suffer from this uncertainty. 

4.1. Transitional Disks 

What is the relevance of MRI accretion to transitional 
disks , i.e., disks with inner holes (e.g ., jE spaillat et al 
20101 [Hughes et ari[2007l [Calvet et al .1[2roE 



35|)'.'' I Chiang fc: 



^ The rotation measures of magnetized FUV-ionized layers could 
Be large, on the order of 10'' rad/m^ at a ~ 30 AU for a facc-on 
HTsk (assuming^ai^jintangled^^ lOmG and an electron colurim 
density of A^ c ~ 3x 10^° cm~''). P erhaps measurements of Faraday 
rotation at radio wavelengths using polarized background active 
galaxies, or polarized emission from the central star itself, coul d Ee 
used to constrain disk magnetic tields. 



Murray-Clay[ (2007, hereafter CMC) proposed that MRI- 
active surface layers at the rim of the hole — i.e., in the 
rim "wall" oriented perpendicular to the disk midplane — 
could supply the observed accretion rates of transitional 
disks. Following CMC, we compute the accretion rate at 
the rim according to the formula 



trans 



3/2 2 
n • 
rim 



3Mrim ^„ fkT\ 

^ 12tt1j X max a x — ^, , 

tdis V M / GM 

(26) 

where Mrim ~ 47ra/iE* {h is the disk half-thickness), 
idiff ~ ^rim/^ time for material to diffuse from 

arim to arim/2, and E* is reinterpreted for transitional 
disk rims as the radial, not vertical, column of MRI- 
active material. We re-compute E* for transitional disk 
rims by taking the grazing angle ^ ^ 1, as is appropriate 
for radiation which penetrates the rim wall at normal 
incidence. 



Figure 10 shows accretion rates for transitional disk 
rims computed acc ording to ( 26 ) , over l aid with obser- 
vations taken from Calvet et aIT f2005). 'Espaillat et al 
([2007| [20081 [20T0| , and [Kim et al. (2009). We have cho- 



sen these references and not others 'because they utilize 
disk models similar enough to each other to yield r eason- 
ably consistent hole radii (cf. 



Merm et al. 2010). Rim 



accretion by FUV ionization seems capable of reproduc- 
ing the trend of increasing Mtrans with increasing flrim- 
However, there remains the problem of transporting the 
material that is dislodged from the rim over the decades 
in disk radius between the rim and the host star. In- 
side the rim, a conventional disk geometry applies, and 
we have already noted in that context that surface layer 
accretion rates at small radii are too small compared to 
those observed. For example, according to Figure [TOj an 
accretion rate of Mtrans ^ lO^^Afo/yr may be initiated 
at Oi-im — 30 AU, but according to Figure [8| this same 
accretion rate cannot be sustained by the MRI inside a 
radius of ^3 AU. Multiple planets co uld solve this pro b- 



lem by shuttling gas inward (PBCll; [Zliu et al.[[2011[ ). 
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